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We have investigated evolution of non-thermal emission from relativistic 
electrons accelerated at around the shock fronts during merger of clusters 
of galaxies. We estimate synchrotron radio emission and inverse Compton 
To . scattering of cosmic microwave background photons from extreme ultraviolet 

(EUV) to hard X-ray range. The hard X-ray emission is most luminous in 



the later stage of merger. Both hard X-ray and radio emissions are luminous 
only while signatures of merging events are clearly seen in thermal intracluster 
medium (ICM). On the other hand, EUV radiation is still luminous after the 
system has relaxed. Propagation of shock waves and bulk-flow motion of ICM 
play crucial roles to extend radio halos. In the contracting phase, radio halos 
are located at the hot region of ICM, or between two substructures. In the 
expanding phase, on the other hand, radio halos are located between two ICM 
hot regions and shows rather diffuse distribution. 



Subject headings: galaxies: clusters: general — hydrodynamics — intergalactic 
medium — particle acceleration — radiation mechanisms: nonthermal 
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1. INTRODUCTION 

Some clusters of galaxies have diffuse non-thermal synchrotron radio halos, which 
extend in a ~ Mpc scale (e.g., Giovannini et al. f993; Rottgering et al. f997; Deiss et al. 
1997). This indicates that there exists a relativistic electron population with energy of a 
few GeVs (if we assume the magnetic field strength is an order of /xG) in intracluster space 
in addition to the thermal intracluster medium (ICM). Furthermore, it is well known that 
such clusters of galaxies have evidences of recent major merger in X-ray observations (e.g., 
Henriksen & Markevitch 1996; Honda et al. 1996; Markevitch, Sarazin, & Vikhlinin 1999; 
Watanabe et al. 1999). In such clusters of galaxies with radio halos, non-thermal X-ray 
radiation due to inverse Compton (IC) scattering of cosmic microwave background (CMB) 
photons by the same electron population is expected (Rephaeli 1979). Indeed, non-thermal 
X-ray radiation was recently detected in a few rich clusters (e.g., Fusco-Femiano et al. 
1999; Rephaeh, Gruber, & Blanco 1999; Kaastra et al. 1999) and several galaxy groups 
(Fukazawa 1999) although their origins are still controversial. In addition to such relatively 
high energy non-thermal emission, diffuse extreme ultraviolet (EUV) emission is detected 
from a number of clusters of galaxies (Lieu et al. 1996; Mittaz, Lieu, & Lockman 1998; 
Lieu, Bonamente, & Mittaz 1999). Although their origins are also unclear, one hypothesis 
is that EUV emission is due to IC emission of CMB. If the hypothesis is right, this indicates 
existence of relativistic electrons with energy of several handled MeVs in intracluster space. 

The origin of such relativistic electrons is still unclear. Certainly, there are point 
sources like radio galaxies in clusters of galaxies, which produces such a electron population. 
However, the electrons can spread only in a few kpc scale by diffusion during the IC cooling 
time. Clearly, this cannot explain typical spatial size of radio halos. One possible solution 
of this problem is a secondary electron model, where the electrons are produced through 
decay of charged pions induced by the interaction between relativistic protons and thermal 
protons in ICM (Dennison 1980). In the secondary electron model, however, too much 
gamma-ray emission is produced through neutral pion decay to fit the Coma cluster results 
(Blasi & Colafrancesco 1999). Moreover, the secondary electron model cannot explicitly 
explain the association between merger and radio halos. 

From N-body -|- hydrodynamical simulations, it is expected that there exist shock 
waves and strong bulk-flow motion in ICM during merger (e.g, Schindler & Miiller 1993; 
Ishizaka & Mineshige 1996; Takizawa 1999, 2000). This suggests that relativistic electrons 
are produced at around the shock fronts through 1st order Fermi acceleration and that 
propagation of the shock waves and bulk-flow of ICM are responsible for extension of radio 
halos. Obviously, the merger shock acceleration model can explicitly explain the association 
between merger and radio halos. However, such hydrodynamical effects on time evolution 
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and spatial distribution of relativistic electrons during merger are not properly considered 
in previous studies. 

In this paper, we investigate the evolution of a relativistic electron population and 
non-thermal emission in the framework of the merger shock acceleration model. We perform 
N-body + hydrodynamical simulations, explicitly considering the evolution of a relativistic 
electron population produced at around the shock fronts. 

The rest of this paper is organized as follows. In §^ we show order estimation about 
the spatial size of radio halos. In ^ we describe the adopted numerical methods and initial 
conditions for our simulations. In §^ we present the results. In §^ we summarize the results 
and discuss their implications. 



2. ORDER ESTIMATION 



In this section, we estimate several kinds of spatial length scale relevant to the extent 
of cluster radio halos. 



2.1. Diffusion Length 

According to Bohm diffusion approximation, a diffusion coefficient is, 

rjEeC 



3eB ' 

where, rj is an enhanced factor from Bohm diffusion limit, is the total energy of an 
electron, c is the velocity of light, e is the electron charge, and B is the magnetic field 
strength of intracluster space. Since IC scattering of CMB photons is the dominant cooling 
process for electrons with energy of ~ GeVs in typical intracluster conditions (Sarazin 
1999), electron cooling time in this energy range is, 

i,, = i.ixioV(|v)" 

where we assume that the cluster redsift is much less than unity. Thus, the diffusion length 
within the cooling time is 

/ 17 \ 1/2 / R \ -1/2 

i^difr~v^~l.lxlO-^Mpc(^^j i^—j . (3) 

This is much less than the spatial size of radio halos. From this reason, electrons which are 
leaking from point sources such like AGNs cannot be responsible for the radio halos. 
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2.2. Shock Wave Propagation Length 



From N-body + hydrodynamical simulations of cluster mergers, the propagation speed 
of shock wave is an order of ~ 1000 km s~^ (Takizawa 1999). Thus, the length scale where 
the shock front propagates during the cooling time of equation is. 



where f shock is the propagation speed of the shock front. Roughly speaking, Lshock.pmp is 
related to the extent of radio halos along to the collision axis since the accelerated electrons 
can emit synchrotron radio radiation only during the cooling time behind the shock front. 



In cluster merger, the shock front spread over in a cluster scale (~ Mpc). Thus, even if 
the shock is nearly standing, radio halos spread in a Mpc scale. Roughly speaking, this is 
related to the extent of radio halos perpendicular to the collision axis. 

In the merger shock acceleration model, therefore, propagation of shock fronts and 
spatial size of shock surfaces play more crucial role to the extent of radio halos than 
diffusion. Furthermore, the model naturally produce 3-dimensionally extended radio halos 
in a Mpc scale. 



We consider the merger of two equal mass (0.5 x IO^^Mq) subclusters. In order to 
calculate the evolution of ICM, we use the smoothed-particle hydrodynamics (SPH) method. 
Each subcluster is represented by 5000 N-body particles and 5000 SPH particles. The initial 
conditions for ICM and N-body components are the same as those of Run A in Takizawa 
(1999). The numerical methods and initial conditions for N-body and hydrodynamical 
parts are fully described in §3 of Takizawa (1999). Our code is fully 3-dimensional. 

To follow the evolution of a relativistic electron population, we should solve the 
diffusion-loss equation (see Longair 1994) for each SPH particle. Since the diffusion term is 
negligible as mentioned in the equation is. 




(4) 



2.3. Spatial Size of Shock Surfaces 



3. 



MODELS 



dN{E,,t) 
dt 



BE, 



d 



e 



[6(Ee,t)iV(Ee,t)]+g(Ee,t), 



(5) 
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where N{Ef,, t)dEe is the total number of relativistic electrons per a SPH particle with 
kinetic energies in the range Ef, to E^, + dE^ (hereafter, we denote kinetic energy of an 
electron to E^,), b{Ee,t) is the rate of energy loss for a single electron with an energy of E^, 
and Q{Ee, t)dEe gives the rate of production of new relativistic electrons per a SPH particle. 

According to the standard theory of 1st order Fermi acceleration, we assume that 
Q{Ec,t) oc E~", where a is described as (r + 2)/(r — 1) using the compression ratio 
of the shock front, r. For the shocks appeared in this simulation, the ratio is roughly 
-\/T0 (Takizawa 1999), which provides a = 2.4. Since it is very difficult to monitor the 
compression ratio at the shock front for each SPH particle per each time step, we neglect 
the time dependence of a. The influence of the changes in a on the results is discussed in 

The normalization of Q{Ee,t) is proportional to the artificial viscous heating, which is 
nearly equal to the shock heating. We generate the relativistic electrons everywhere even if 
explicit shock structures do not appear in the simulation. We assume that sub-shock exists 
where a fiuid element has enough viscous heating. Such sub-shocks are recognized in higher 
resolution simulations (e.g. Roettiger, Burns, & Stone 1999). We assume that total kinetic 
energy of accelerated electrons from ii^c = to +oo is 5% of the viscous energy, which is 
consistent with the recent TeV gamma ray observational results for the galactic supernova 
remnant SN 1006 (Tanimori et al. 1998; Naito et al. 1999). Note that equation ^ for the 
evolution of a relativistic electron population is linear in N{Ee, t). Thus, it is easy to rescale 
our results of N{Ei,,t) if we choose other parameters for the acceleration efficiency. We 
neglect energy loss of thermal ICM due to the acceleration. 

For b{Ee,t), we consider IC scattering of the CMB photons, synchrotron losses, 
and Coulomb losses. We neglect bremsstrahlung losses for simplicity, which is a good 
approximation in typical intracluster conditions (Sarazin 1999). Then, if we ignore weak 
energy dependence of coulomb losses, the loss function b{Ee,t) becomes 



where bc{t) = 7.0 x 10-^%n^{t)/cm-^) and bi = 2.7 x 10^^^ + 2.6 x 10-^^{B{t)/ fiGf 
(if b{Ef,,t) and E^ are given in units of GeV s^^ and GeV, respectively.). In the above 
expressions, Uf, is the number density of ICM electrons and B is strength of the magnetic 
field. 

To integrate equation (|^) with the Courant and viscous timestep control (see Monaghan 
1992), we use the analytic solution as follows. First, we integrate equation (^ from t to 
t + At, regarding the second term on the right-hand side as being negligible small. Then, 



b{E,,t) = b,{t)+bi{t)El 



(6) 




(7) 
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where 



Eefl = Y — tan|^ arctan J j^^e + ybcbiAt 

Next, we add the contribution from the second term to the above N{Ec,t + At) using the 
second-order Runge-Kutta method. In the present simulation, N{Ee,t) is calculated on 
logarithmically equally spaced 300 points in the range Eg = 0.05 to 50 GeVs for each SPH 
particle. 

Magnetic field evolution is included by means of the following method. We assume 
initial magnetic pressure is 0.1 % of ICM thermal pressure. This corresponds to -B = O.l/xG 
in volume- averaged magnetic field strength. For Lagrangean evolution of B, due to the 
frozen-in assumption we apply B(t)/B{tQ) = (picm(^)/picm(^o))^''^- Field changes due to 
the passage of the shock waves is not considered in this paper. The change may depend on 
field configuration at the shock front and have value of ~ 1 — 4. However, it is difficult to 
examine it in the present simulation even under high f3 condition. We will try this problem 
in the future paper. 

Our model implies continuous production of power law distributed relativistic electrons 
at around the shock fronts. This is valid only when Atacc is sufficiently shorter than the 
dynamical timescale of the system (~ 10^ yr), where Atacc denotes acceleration time in 
which Q{Ee,t) becomes power law distribution. It is presented in the framework of the 
standard shock acceleration theory as Atacc = Sru~^{r — + rK2) , where u is the 

flow velocity of the upstream of the shock front, and /ti^2 are diffusion coefficients of the 
upstream and downstream, respectively (see e.g. Drury 1983). Assuming Bi = B2 and 
Bohm diffusion approximation of equation (|l|), 

^ VGeVy V102AlO%m s^V \fiGj ^ ' 

This value is certainly much shorter than the dynamical timescale. 



4. RESULTS 

Figure |1] shows the time evolution of non-thermal emission for various energy band: 
from top to bottom, IC emission of the Extreme Ultraviolet Explorer (EUVE) band (65-245 
eV), soft X-ray band (4-10 keV), and hard X-ray band (10-100 keV), and synchrotron 
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radio emission (10 MHz - 10 GHz). The times are relative to the most contracting epoch. 
The calculation of the luminosity for each band is performed in the simplified assumption 
that electrons radiate at a monochromatic energy given by Ex = 2.5keV(£'c/GeV)^ and 
u = 3.7MHz(i?//iG)(i?e/GeV)^ for IC scattering and synchrotron emission, respectively. 
Since cooling time is roughly proportional to in these energy range, the higher the 
radiation energy of IC emission is, the shorter duration of luminosity increase is. In other 
words, luminosity maximum comes later for lower energy band. Hard X-ray and radio 
emissions come from the electrons with almost the same energy range. The luminosity 
maximum in the hard X-ray band, however, comes slightly after the most contracting epoch. 
On the other hand, radio emission becomes maximum at most contracting epoch since the 
change of magnetic field due to the compression and expression plays an more crucial role 
than the increase of relativistic electrons. In any cases, radio halos and hard X-ray are 
well associated to merger phenomena. They are observable only when thermal ICM have 
definite signatures of mergers such as complex temperature structures, non-spherical and 
elongated morphology, or substructures. Soft X-ray emission, which is observable only in 
clusters (or groups) with relatively low temperature (~ 1 keV) ICM, is still luminous in ~ 
1 Gyr after the merger. Thus, the association of mergers in this band is weaker than in the 
hard X-ray band. Moreover, EUV emission continues to be luminous after the signatures of 
the merger have been disappeared in the thermal ICM. 

Figure |^ shows the IC spectra at t = 0.0 (solid lines) and 0.25 (dotted lines). In lower 
energies, L^, oc z/~°-^, which is originated from the electron source spectrum Q{Ee) oc -E"^'^- 
On the other hand, in higher energies, the spectrums become close to steady solution, 
Lj, oc z/~^'^, owing to the IC and synchrotron losses (Longair 1994). The break point of the 
spectrum moves toward lower energies as time proceeds. 

Figure |^ shows the snapshots of synchrotron radio (lOMHz-lOGHz) surface brightness 
distribution (solid contours) and X-ray one of thermal ICM (dashed contours) seen from the 
direction perpendicular to the collision axis. Contours are equally spaced on a logarithmic 
scale and separated by a factor of 7.4 and 20.1 for radio and X-ray maps, respectively. 
At t = —0.25, the main shocks are located between the two X-ray peaks and relativistic 
electrons are abundant there. Thus, the radio emission peak is located between the two 
X-ray peaks although the magnetic field strength there is weaker. At t = 0.0, relativistic 
electrons are concentrated around the central region since the main shocks are nearly 
standing and located near X ~ ±0.2. Furthermore, gas infall compress ICM and the 
magnetic field. Thus, radio distribution shows rather strong concentration. In these phase 
(at t = —0.25 and 0.0), the radio halo is located at the high temperature region of ICM. On 
the other hand, at t = 0.25, relativistic electron distribution becomes rather diffuse since 
fresh relativistic electrons are producing in the outer regions as the shock waves propagate 
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outwards. At t = 0.25 the main shocks are located at X ~ ±1. Between the shock fronts 
rather diffuse radio emission is seen. In this phase, the radio halo is located between two 
high temperature regions of ICM. 

As described above, the morphology of the radio halo is strongly depending on the 
phase of the merger when viewed from the direction perpendicular to the collision axis. 
When viewed nearly along the collision axis, however, this is not the case. Figure ^ shows 
the same as figure |^, but for seen from the direction tilted at an angle of 30° with respect 
to the collision axis. Radio and X-ray morphology are similar each other in all phases. 
When the cluster is viewed along the collision axis, the distribution of relativistic electrons 
roughly follows that of the thermal ICM since the shock fronts face to the observers and 
spread over the cluster. The distribution of magnetic field strength also roughly follows 
that of the thermal ICM. Therefor, the radio morphology follows X-ray one. 

Figure ^ shows the synchrotron radiation spectra at t = 0.0 (solid lines) and 0.25 
(dotted lines). At t = 0.0, the synchrotron spectrum follows that of IC emission. On 
the other hand, at t = 0.25 there exits a bump at lower energies in the spectrum, which 
cannot be seen in that of IC emission. Emissivity of synchrotron radiation depends on not 
only relativistic electron density but also magnetic energy density, which is larger in the 
central region in this model. On the other hand, IC emissivity depends on the electron 
density and CMB energy density, which is homogeneous. Thus, the total synchrotron 
spectrum is more like that in the central region than the total IC spectrum. At t = 0.0, 
since relativistic electrons are centered, the emission from the outer region is negligible for 
both spectra. Thus, similar results are given. On the other hand, at t = 0.25, propagation 
of shocks makes a diffuse distribution of relativistic electrons. Therefor, the contribution 
from the outer region is not negligible for the total IC spectrum while the total synchrotron 
spectrum is still biased the central region as seen in figure ^. In the central region, however, 
electrons produced by the main shocks at t ^ with energies more than ~ several GeVs 
have already cooled down. Thus the spectrum in the central region have a bump in lower 
energies, which is present in the total spectrum. The emission above ~ 30 MHz is mainly 
due to the electrons produced by the sub-shocks there. If such sub-shocks do not exist, the 
emission near the main shocks, where the magnetic field strength is rather weaker, can be 
seen in this energy range. Note that this feature of the synchrotron spectrum is sensitive to 
the spatial distribution of magnetic field. 
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5. CONCLUSIONS AND DISCUSSION 

We have investigated evolution of non-thermal emission from relativistic electrons 
accelerated at around the shock fronts during merger of clusters of galaxies. Hard X-ray and 
radio radiations are luminous only while merger signatures are left in thermal ICM. Hard 
X-ray radiation becomes maximum in the later stage of merger. In our simulation, radio 
emission is the most luminous at the most contracting epoch. This is due to the magnetic 
field amplification by compression. According to the recent magnetohydrodynamical 
simulations (Roettiger, Stone, & Burns 1999), however, it is possible that the field 
amplification occurs as the bulk flow is replaced by turbulent motion in the later stages of 
merger. If this is effective in real clusters, radio emission can increase by a factor of two or 
three than our results in the later stages of merger. EUV emission is still luminous after 
the merger signatures have been disappeared in thermal ICM. This is consistent with the 
EUVE results. 

Morphological relation between radio halos and ICM hot regions is described as follows. 
In the contracting phase, radio halos are located at the hot regions of thermal ICM, or 
between two substructures (see the left panel of figure |^). This may correspond to A2256 
(Rottgering et al. 1994). In the expanding phase, on the other hand, radio halos are located 
between the two hot regions of ICM and show rather diffuse distribution (see the right panel 
of figure This may correspond to Coma (Giovannini et al. 1993; Deiss et al. 1997) and 
A2319 (Feretti, Giovannini, & Bohringer 1997). In the further later phase, the shock fronts 
reach outer regions and the GeV electrons are already cooled in the central parts. Then, 
radio halos are located in the cluster outer regions near the shock fronts and we cannot 
detect radio emission in the central part of the cluster. However, observational correlation 
between ICM hot regions and radio halos is not clear since the electron temperature there is 
significantly lower than the plasma mean temperature due to the relatively long relaxation 
time between ions and electrons (Takizawa 1999, 2000). Note that until now we could only 
find the electron temperature through X-ray observations. This may correspond to A3667 
(Rottgering et al. 1997). It is possible that such radio 'halos' located in the outer regions 
are classified into radio 'relics' since their radio powers and spatial scales becomes weaker 
and smaller than those of typical radio halos, respectively. When the cluster is viewed 
nearly along the collision axis, however, such morphological relations between radio halo 
and ICM are unclear and radio and X-ray distributions become similar each other. 

We neglect the changes of the spectral index in the electron source term. Since the 
mach number is gradually increasing as merger proceeds (Takizawa 1999), the spectrum of 
relativistic electrons becomes flatter as time proceeds. We believe that such changes in the 
spectral index does not influence our results seriously because most of relativistic electrons 
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are produced in the central high density region, where the mach number is almost constant. 
In the later stage of merger, however, contribution of relativistic electrons produced in the 
outer region cannot be negligible in higher energy range (~ 10 GeV) since cooling time is 
relatively short. Thus, it is probable that the inverse Compton spectrum in the hard X-ray 
(~ 10 — 100 keV) becomes flatter in the later stages. 

The lower energy part of the electron spectrum can emit hard X-ray through 
bremsstrahlung. Whether IC scattering or bremsstrahlung is dominant in the hard X-ray 
range is depending on the shape of the electron spectrum. Roughly speaking, when the 
spectrum of relativistic electrons is flatter than IC scattering dominates the other and 

vice versa (see Appendix). Furthermore, the shape of the electron spectrum in the lower 
energy part is flatter than the originally injected form since the cooling time due to the 
coulomb loss is proportional to and very short (Sarazin 1999). In the present simulation, 
therefor, it is most likely that the contribution of the bremsstrahlung components in the 
hard X-ray range is negligible. More detailed calculations, including nonlinear effects for 
the shock acceleration (Jones & Ellison 1991), by Sarazin & Kempner (1999) show that 
IC scattering is dominant in the hard X-ray when the accelerated electron momentum 
spectrum is flatter than p^^'"^ , which corresponds to the electron energy spectrum of 
in the fully relativistic range. Thus, the bremsstrahlung contribution in the hard X-ray 
should be considered in mergers with low mach numbers. 

A merger shock acceleration model also predicts some gamma-ray emission. Electrons 
which radiate EUV due to IC scattering also emit ~ 100 MeV gamma-ray through 
bremsstrahlung. Furthermore, it is most likely that protons as well as electrons are 
accelerated at around the shock fronts. Such high energy protons also produce gamma-rays 
peaked at ~ 100 MeV through decay of neutral pions. Although the energy density ratio 
between electrons and protons in acceleration site is uncertain, the contributions of protons 
and the bremsstrahlung to the emission become important in the higher energies. We think 
that it is interesting to investigate from hundreds MeV to multi TeV emissions, which 
are observable with operating instruments such as EGRET, ground-based air Cerenkov 
telescopes, and planning projects like GLAST satellites. However, since the diffusion length 
of protons within the cooling time is much longer than that of electrons, treatment of the 
diffusion-loss equation for protons are more complex than the model in this paper. 

We would like to thank Drs. Y. Fukazawa and S. Shibata for helpful comments. MT 
thanks Drs. S. Mineshige and T. Shigeyama for continuous encouragement. MT is also 
grateful to S. Tsubaki for fruitful discussion. 
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A. ESTIMATION OF THE BREMSSTRAHLUNG CONTRIBUTION IN 

THE HARD X-RAY RANGE 

Wc estimate the contribution of the bremsstrahhmg form the lower energy part of 
the electron spectrum to the hard X-ray emission in our model, which is neglected in this 
paper. Although the crude estimations discussed here are order-of magnitudes, it is helpful 
to explain which mechanism, IC scattering and bremsstrahlung radiation, is dominant. 

The standard theory of 1st order Fermi acceleration provides power law spectrum in 
momentum distribution. We, therefore, assume the momentum spectrum of accelerated 
electrons has a form of 

= cNo ( ) [cm-3 (eV 0-^'] , (Al) 



dPe XlTbcC 

where is an electron momentum, c is the speed of light, and rUe is the electron 
rest mass. Using the relation between the momentum Pe and the kinetic energy E^, 
Ee = {P^c^ + mlc^Y^"^ — rrieC^, the spectrum for kinetic energy is given by 



dEe nieC^ KirieC 

For relativistic electrons, Ee ^ mgC^, the spectrum becomes 



^^==iV„fA) . (A3) 



dEf, KUleC 



which is consistent with an assumption of Q{Ee,t) in section 3. Since the hard X-ray 
photons are produced from CMB photon field via IC scattering of relativistic electrons 
in our model, we use this form for the estimation of IC emissivity. For non-relativistic 
electrons, E^ <^ rrieC^, the spectrum becomes 



^^2-^N„?l±^(^y''' , (A4) 

where we use the relation Ee = Pg^/(2me). Since the bremsstrahlung radiation at the hard 
X-ray range is emanated from electrons with almost the same energy range, we use this 
form for the estimation of bremsstrahlung emissivity. 

For simplicity, we approximate the emissivity, e, for both IC and bremsstrahlung 
processes to be, 

AAT rIT? rITP 

erg s"^ cm"^ eV""^] , (A5) 



dNe dEe 


dEe 


dEe ds^ 


dt 
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where the emission rate is assumed to be equal to the electron energy loss rate and is the 
photon energy. 

To obtain the bremsstrahlung emissivity, we assume that an electron with energy Eg 
loses its energy to emit photon of energy = Ee after it has traversed one mean free path 
Xq. Hence, 



dE, 

de^ 



(A6) 



and 



dt 



Vr, 



erg s" 



(A7) 



where is electron velocity. From aTn^XQ ~ 1, we approximate 



dt 



where ctt denotes the cross section of Thomson scattering and no denotes the density 
of ambient matter. Using equations (|A4|) , (|A5|), ( |A6| ) , and ([A8|) , we estimate the 
bremsstrahlung emissivity in the hard X-ray energy range, = £:hxr; as 



Cbrem ~ 2 2 NoaTCnQTrieC 1 H 



\ / \ — —+1 

^HXR \ / ^HXR \ 2 ^ 



(A9) 



For IC scattering of CMB photons, the photon energy after the scattering by an 
electron with energy E^ = me7eC^ is approximated by single energy of = 7e^£^cMB, where 
^CMB is peak energy of CMB spectrum. Thus, 



dE^ _ 1 /£cmb\ 5 / e 
de^ 2 V mpC^ / v 



(AlO) 



Since the scattering is in Thomson energy range (ecMB7e ^ meC^), we set 



dE, 



dt 



-crTC7e^ecMB^CMB [erg s ^] 



(All) 



where ^cmb is photon number density of CMB fields. Using equations (|A^), ([A^), ( [A10| ), 
and ( All ), we estimate the IC emissivity in the hard X-ray energy range, = £:hxr; as 



O / — \ , ot — \ 

2,^ _ /£^CMB\ 2 /enxR \ 2 

eiC — gJVoCTTCecMB^CMB 



(A12) 



- 13 - 



From equations (^) and ( |A12|) , we derive the emissivity ratio as 



Cbrem _ 3 Uq / ^ £hXr\ / £cMb\ 2 / ^HXr X 2 /^ gCMB \ 



eic 2ncMB V rripC 
For ebrcm < eiC) we get the relation 



mpC^ 



2 In 



a < 



2 ncMB 



]^ _|_ ^hxr A (£cmb£hxr) 



moC- 



ln( 



2ecMB ] 



(A13) 



(AM) 



Substituting typical values for ICM at z ~ 0, = 1 x 10^^ cm~^, wcmb = 400 cm~^, 
£cMB = 6.57 X 10"^ eV, and ^hxr = 10 keV, we obtain a < 2.5. At ^hxr = 100 keV, the 
spectrum of equation ( |A9| ) becomes steeper as electrons enter the trans-relativistic energy 
range. As a result, the bremsstrahlung emissivity is reduced so that the larger index of 
electron spectrum is accepted for ebrem < cic- 
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(c) 10-100 keV 
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(d) 10MHz -10GHz 
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Fig. 1. — The time evolution of non-thermal emission for various energy band: from top to 
bottom, inverse Compton scattering of EUVE band (65-245 eV), soft X-ray band (4-10 keV), 
and hard X-ray band (10-100 keV), and synchrotron radio emission (10 MHz - 10 GHz). The 
times are relative to the most contracting epoch. 
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Fig. 2. — The inverse Compton scattering spectra at t — 0.0 (solid lines) and 0.25 (dotted 
lines) . 
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Fig. 3. — Snapshots of synchrotron radio (lOMHz-lOGHz) surface brightness distribution 
(sohd contours) and X-ray one of thermal ICM (dashed contours) seen from the direction 
perpendicular to the collision axis. Contours are equally spaced on a logarithmic scale and 
separated by a factor of 7.4 and 20.1 for radio and X-ray maps, respectively. 




Fig. 4. — Same as figure ^, but for seen from the direction tilted at an angle of 30° with 
respect to the collision axis. 
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Fig. 5. — The synchrotron radiation spectra at t = 0.0 (soUd hnes) and 0.25 (dotted hnes). 



